Additivity properties of a Gaussian channel 
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The Amosov-Holevo- Werner conjecture implies the additivity of the minimum Renyi entropies at 
the output of a channel. The conjecture is proved true for all Renyi entropies of integer order greater 
than two in a class of Gaussian Bosonic channel where the input signal is randomly displaced or 
where it is coupled linearly to an external environment. 



PACS numbers: 03.67.Hk,03.67.-a,03.65.Db,42.50.-p 



o 
o 

c 
3 



> 
in 

o 
m 
o 

o 

^— > ■ 
G ■ 
cd 

3 : 
cr 



One of the most challenging open questions of quan- 
tum communication theory is the additivity of the var- 
ious quantities characterizing the information transmis- 
sion in a channel Q] . The issue at hand is whether quan- 
tum entanglement is able to improve the performance 
of classical protocols @>,0- The supposed additivity of 
the Holevo information |4j is the most important exam- 
ple of this kind of issue. The maximum of this quan- 
tity over all possible encoding procedures is known to 
provide the capacity C\ in transmitting classical infor- 
mation for a single use of the channel. However, if the 
sender of the message is allowed to encode messages in 
entangled states among m successive uses of the commu- 
nication line, then the resulting capacity per channel use 
might be higher than C\ Q. For this reason to com- 
pute the ultimate classical capacity C of the channel it 
is necessary to introduce a regularization of the Holevo 
information where a limit m — > oo has to be performed 
P, HJ . All this could be avoided if only the Holevo infor- 
mation was shown to be an additive quantity. Up to now 
no channel has been found for which this regularization 
is necessary: on the contrary all the channels for which 
the value of C has been calculated have additive Holevo 
information HHH)]. 

The additivity of the Holevo information has been 
linked to the additivity of other relevant quantities in 
Ref. Q. In particular, it is known that proving the addi- 
tivity of the Holevo information is equivalent to proving 
the additivity of the minimum von Neumann entropy § 
at the output of the channel. Given a channel described 
by the completely positive (CP) linear map M. on the 
input space TL, this quantity is defined as 
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S(M® m {p)) 



(1) 



where the minimization is performed over all the possi- 
ble input states p of m successive uses of the channel, 
and where S(p) = — Tr[plnp]. The additivity hypothesis 
requires §(A4® m ) to be equal to m times the minimum 
entropy for a single use of the channel S(A4): this con- 
jecture seems simpler to study than the additivity of the 
Holevo information and some authors have focused their 
attention to it [a, l2l llOL 111) . As a matter of fact the 
alleged additivity of the S is just a particular instance 
of the Amosov-Holevo- Werner conjecture ^3 which re- 



quires the maximum of the output z-norm v z (M® m ) of 
the channel to be multiplicative, i.e. it requires that for 
m integer 



v z {M 



(g>m\ 



max \\M® m (p)\\ z = [v z {M)] m (2) 
pew®" 1 



where the maximization is performed again over all the 
input states of m uses of the channel and where 



\\A\\ z = (Tr\An 1/z Ol 



(3) 



is the z-norm of the operator A. In other words, the 
conjecture requires the maximization in the left-hand- 
side of Eq. (J2J to be achieved on non entangled states of 
H® m . The connection between the property of Eq. (J2J 
and the additivity of the minimal output entropy, can be 
established through the quantum Renyi entropy 



S z (p) 



\nTr[p z 
z-1 



(4) 



This quantity is monotonic with respect to the z-norm 
© of the state p. For z — 2 the Renyi entropy is a 
function of the linear entropy 1 — Tr[/9 ] and in the limit 
z — > 1 it tends to the von Neumann entropy [lflj . As for 
the case of S one can define the minimal value 
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S z (M® m (p)) 



(5) 



If the Amosov-Holevo- Werner conjecture (J2J is true then 
the minimum output z-Renyi entropy is additive and vice 
versa. Moreover, if such property is verified for values of 
z arbitrarily close to 1 then the additivity of § (and hence 
of the Holevo information) follows p^j . 

In this paper we will analyze the conjecture @ for a 
set Gaussian channels and prove that it is true for all 
integer z. The material is organized as follows. In Sect. 
[I] we introduce the simple Gaussian channel model Af n 
and in Sees. II Al and ll Bl we show that the conjecture (j2J) 
applies to this channel when z is integer. In Sec. II CI 
we analyze the case of generic z giving some bounds for 
v z {N® m }. In Sec. HI1 we generalize the results of the first 
section to a whole class of Gaussian channels. 
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I. THE CHANNEL MODEL 



A. The proof 



The channel we analyze here is a Bosonic linear channel 
where the photonic signal from the sender is displaced 
randomly by the environment. This system is described 
by the CP map Af n which transforms the input state of 
the channel into the output 



Af n (p) = / d 2 /i P n (/i) D<ji)pD^{n) 



(6) 



where, for n 0, P n (fi>) is the circularly symmetric prob- 
ability distribution 



(7) 



and D(p) = exp{p,a* — P*a) is the displacement operator 
of the annihilation a of the input signal. This channel is 
Gaussian, i.e. it maps the set of input states with Gaus- 
sian symmetrically characteristic function into itself |l5j | . 
Moreover, the map JBJ is unital (i.e. it transforms the 
identity operator in to itself) and it is covariant under 
displacement or phase transformation [ljj. When Af n 
acts on a coherent state p a = \a)(a\ the following trans- 
formation takes place, 



with 



p a -> N'niPa) = D(a) r(n) D\a) 



r(n) 



1 



1 



(8) 



(9) 



the thermal state that gives the output of the channel for 
a vacuum input [ljj. The state N n {p a ) has z-norm © 
equal to 



\Wn{Pa)\\ 



1 



(n + l) z — n z 



1/z 



(10) 



which does not depend on a since it is invariant under 
the unitary transformation D(a). In Ref. | ll| the right- 
hand-side of Eq. (fTTTf) was shown to coincide with the 
z-norm of the single use of the channel v z {M n ): at least 
for all z = k integer. In Sec. II Al we will generalize this 
result showing that, for all integer k, the classical channel 
satisfies the identity, 



1 



- infk 



{n + iy 



(ii) 



hence proving the conjecture J3J) for integer z = k for 
the channel Af n . Equations (|1(J|> and (|llf> imply that the 
maximization implicit in the definition of ^(A/ - ® m ) is 
achievable with separable input states of the form \ai) i <g> 
•••(g) \a m ) m , i.e. by feeding the channel with a coherent 
state in each of the m successive uses. This result will be 
proved explicitly in Sec II Bl 



In this section we show that Eq. (|1 1|) applies for inte- 
ger z. Clearly, the right-hand-side of this equation is a 
lower bound for the left-hand-side: the former is in fact 
the output z-norm associated to the input signal where 
the m uses of the channel have been prepared in coher- 
ent states. To prove the equality in Eq. I|llfl it is hence 
sufficient to show that the right-hand-side is also an up- 
per bound for ^(A"®" 1 ), i.e. that for all input states 
p G 7^® m the following inequality applies, 



Tr{[A;® m (p)] fc }< 



1 



{n + iy 



(12) 



The method to derive this property is similar to the one 
given in Ref. ^l| where an analogous approach was used 
to calculate the minimum output Renyi entropy I0J of 
integer order for a single channel use (m =1). The only 
difference is that here we are dealing with an extra ten- 
sorial structure associated with m > 1. For the sake of 
clarity we divide the proof in two separate parts. First we 
show that the quantity on the left-hand-side of Eq. l(T2^ 
can be expressed as the expectation value of an diago- 
nalizable [23] operator which acts on the Hilbert space 
(7i® m )® fc : this allows us to derive an upper bound for 
Tr {[N~n m (p)] k } by considering the eigenvalue A of 6 
with maximum absolute value. The second part of the 
proof is devoted to the analysis of the tensorial struc- 
ture of and to the proof that Ao coincides with the 
left-hand-side of Eq. ljT2|) . 

Part one:- Without loss of generality we can assume 
the initial state of the m uses of the channel to be pure, 
i.e. p = \il)}(ip\. The convexity of the norm @ guaranties 
in fact that the maximization in Eq. (0) is achievable 
with pure input states 0, : thus if Eq. (|T2"|) holds 
for all pure states, then it is valid also for all the other 
channel inputs. In general will be entangled among 
the various channel uses and the corresponding output 
state will be 

A^f m (p) = J dVl ' ' ' d 2 MnA(Ml) ' ' ' Pn(Prn) 

DM ■ ■ ■ D m (fx m ) p D\p x ) ■ ■ ■ Dl{p m ) , (13) 

where D r {p) = exp(pal — p*a r ) is the displacement as- 
sociated with the annihilation operator a r of the rth use 
of the channel. Equation (|13fl can be expressed in a 
more compact form by introducing a vectorial notation, 
where p, = (pi, • • • , p m ) is a complex vector in C m and 
a = (ai, • • • , a m ). The output state becomes thus 

A^? m (p) = / d 2 p P n (p) D(p) p D\p) , (14) 



where 



exp[-|/2| 2 /«] 



(15) 
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and D(p) = exp(/2 • — a ■ fp) is a multi-mode displace- 
ment operator where the input of the rth use of the chan- 
nel is displaced by p r . Consider now for z — k integer 
the quantity 

Tr{K® m (p)] fc } = / d 2 pi---d 2 p k P n (p 1 )---P n (p k ) 



(16) 



Since p = the trace in the integral can be ex- 

pressed as the expectation value of an operator which 
acts in a extended Hilbert space (7^® m )® fe made of k 
copies of the initial one. In fact, from the invariance of 
the trace under cyclic permutation we have 



TrlDfafrDifa) ■ ■ ■ D{p k )pD\p k )] = (^\D\^)D(p 2 )\i;)(ij\D\^)D(^M ■ ■ ■ (^04 Wi)|V>> 

= Tr {(pOp®---®/)) D\([li)Di(fl2)®Dl(p2)D2{fl3)^---^Dl(ll k )Dk(ixi) } , (17) 

I 



where the k scalar products in the input Hilbert space 
7i® m in the first line were replaced with a single expec- 
tation value in ( / ft® m )® fc m the second line. In Eq. (TTJ 
the operator D s {p) represents the multi-mode displace- 
ment that operates on the sth copy of 7i® m , i.e. 



D s (p) = exp(p-al 



a s ■ p 



where for s = 1, 



■ k 



(18) 



(19) 



are the m annihilation operators pertaining to the sth 
copy of Ji® m . With this trick Eq. (fl?)|) can be written as, 

Tr {\Mf n (p)] k } = ■ • <8> p) 0] , (20) 

where each of the k copies of the state p is associated 
to one of the multi-mode annihilation operator a s and 
where is the operator on (j-(® m ^® k given by 



(21) 



9 = / d pi ■ ■ ■ d p k P n (pi) ■ ■ ■ Pn(pk) 



x D\{pi)Di{p2) ® • • • ® Dlift^Dkipt) . 

Equation (|20|) allows us to derive an upper bound for the 
quantity on the left-hand-side by considering the eigen- 
value Aq of with maximum absolute value, i.e. 



Tr{[A/r>)] fc H|Ao| 



(22) 



Part two:- To calculate Ao it is useful to analyze in 
details the properties of the operator 0. As shown in 
App. IA II this operator has a very simple tensorial form 
with respect to the index r. In fact Eq. I|22|l can be 
decomposed as 



= (g) r , 



(23) 



r=l 



where, for r — 1, ■ ■ • ,m, the operator r acts on the 
modes associated with the annihilation operators 



a r— (air, a2 r , ■ • • , Q kr ) . 



(24) 



In vectorial notation r can be expressed as 
d 2 P, 



Q r = 



(25) 



where, as in Eq. (|24|) . p r = (pi r , P2r, • ■ • , Pkr) is a k- 
element vector, and where G and 



C 



1 

n 



A 
~2 



(26) 



are k x k real matrices (II is the identity). For k ^ 3, G 
and A are 



G = 







A = 



-1 1 
-1 1 
0-1 





1 

0-1 

1 -1 

1 




-10 







-1 1 

-1 

1 



-1 

1 



(27) 



(28) 



For k — 2 A is null, while for k = 1 both G and A are 
null. The decomposition 112MI) shows that in the right- 
hand-side of Eq. (|20|l we have a product of two operators 



of (7i® m )® k which have 
sition: the operator p ® 
the index r = 1, ■ • • , k, 
to the index s = 1, • • • 
sive uses of the channel 



'orthogonal" tensorial decompo- 
■ ■ ® p factorizes with respect to 
while factorizes with respect 
, m associated with the succes- 
This property is common to all 
memoryless channels where the corresponding CP map 
acts on each channel use independently. However, the 
channel model we are considering allows us to further 
decompose the operator 0. In fact, G and A of Eqs. I|27|) 
and (|28|) are two circulant matrices |l6j | which commute 



4 



and possess a common basis of orthogonal eigenvectors. 
This means that there exists a k x k unitary matrix Y 
for which 



D = YCY^ 
E = YGY^ 



1/n + Y AY^ 



(29) 



are diagonal. Since A is antisymmetric, its eigenvalues 
i^j are imaginary and the diagonal elements of D (i.e. 
the quantities dj = l/n + i£j) have positive real part. 
Using these properties we can rewrite the operator Q r 
of Eq. (|25j) in factorized form by performing the change 

of integration variables v r =^ r -Y* and introducing the 
new annihilation operators 



b r = (bir,b 2r , • • • , b kr ) =a r -Y* 
These operations yield 



Qr = 



(30) 



(31) 



with 



Qj r — 



Ip/^e-^f^H, (32) 

where D\, jr {y) = exp[^&j r — u*bj r ] is the displacement 
operator associated with bj r , while ej is the jth diagonal 
elements of the matrix E (i.e. the jth eigenvalue of G). 
As demonstrated in Ref. 0, 0], this expression can 
be further simplified, proving that Qj r is diagonal in the 
Fock basis of the mode bj r and equal to 



2/n 



2di 



2d, 



2d 3 + \e J \ 2 



(33) 



(for the sake of completeness we give an alternative 
derivation of this result in App. IA 2(1 . 

Equations (|23() . 1|31|) and i|33[l show that the eigenval- 
ues of Q are products of eigenvalues of Qj r . In particular 
Ao of Eq. ll'L'l) is obtained by taking the eigenvalues of 
the Qj r that have maximum absolute value. Since the 
constants dj have positive real part, the quantities we 
are looking for are 2/[n(2dj + \ej\ 2 )], i.e. they are the 
eigenvalues of the operators Qj r associated with the vac- 
uum state of the mode bj r . This allows us to express the 
value of Ao as 



An 



nn 

r=l j=l 



2/n 



2d, 



A' 

det[C + GtG/2] J 



{n + l) k 



(34) 



which, replaced in Eq. I)22[l. proves the thesis i|12jl . [In 
deriving the second identity we have used the invariancc 
of the determinant under the unitary transformation Y, 
while the last identity can be obtained from the defini- 
tions <|27|l and l|28|l by direct calculation of the determi- 
nant itself]. 



B. Optimal inputs 

Equations lfTU|) and ((TTj) prove that tensor products of 
coherent states are optimal since they allow the channel 
J\f n to achieve the maximal A;-norm at the output of m 
successive uses. Here we re-derive this result by showing 
that the state p (g> • • • ® p of Eq. (|20|l with p given by a 
tensor product of coherent states in the inputs modes a s , 
is an eigenvector of 9 associated with the eigenvalue Ao 
of Eq. JSJl. 

From the analysis of the previous section we know that 
the eigenvectors of Aq can be written as 



(35) 



where |$ r ) r is an eigenvector of 9 r of Eq. I|31(l relative 
to its eigenvalue with maximum absolute value. For in- 
stance, in deriving Eq. (|3*ljl we have considered the state 

where each of the br modes is in the vacuum. However 
this is not the only possibility. In fact, we notice that for 
any k the matrices G and A of Eqs. 127|) and 128(1 have a 
null eigenvalue (say for j — 1), associated with the com- 
mon normalized eigenvector (1, 1, • • • : l)/^/k. On one 
hand this means that all the elements in the first row 
of the matrix Y of Eq. (|29|l are equal to l/\/k. On the 
other hand, this implies also that e\ = 0, d\ = l/n and, 
according to Eq. (J2SJ), 



Qir — lir , ■ 
This means that any state of the form 

|$ r ) r = \<f> r )b lr ® |0)fc 2r <g> • •• (g> \0) bk 



(36) 



(37) 



where the mode b\ r is prepared in a generic state | <p r ) 
while the other bj r are in the vacuum, is an eigenstate 
of O r relative to the eigenvalue with maximum absolute 
value. Consider now the case of \<j) r ) — \\fka r ) coherent. 
By usin g th e symmetric characteristic function decom- 
position [18j we can show that, when expressed in terms 
of the operators a r , the state (|37|l is a tensor product of 
coherent states \a r ). In fact, defining the complex vector 

7= (Vka r , 0, • • • , 0) we can express the state |$ r ) r as, 



d 2 v 



exp[-| P| 2 /2 



+ V -{b r ~ 7) f - (br - 7)' V 



d 2 V 



exp[-| P | 2 /2 



+ p-(a r -Y^l)-{a r -l-Y)-p\ 



\0h)a lr { a r\ 



r/a k , 



(38) 



where the second identity is obtained by substituting the 
v with P -Y> in the integral, while the third identity de- 
rives from the properties of the matrix Y discussed above. 
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The thesis finally follows by replacing this expression in 
Eq. 



l*X*l = ®(\<*r)*M 
r=l 



K)a fc >r| ) ,(39) 



and observing that this can be represented as p® ■ ■ ■ ® p 
of Eq. CSJ with p= ((g)™ 1 |a r )(a r |)® fe . 



This bound must be compared with the lower bound 



(n + 



m/z 



(43) 



for arbitrary z that derives by considering as input of 
the m successive uses of the channel a tensor product of 
coherent states. This two bounds are plotted in Fig. 



C. Upper bound 



II. GENERALIZATION 



In this section, starting from the values of the v z {N® m ) 
for z integer derived in the previous section, we give some 
upper bounds for the channel z-norm of generic order. 




FIG. 1: Plot of the bounds of the 2-norm v z (N® m ) of the 
channel Af n . The function v z (J\f® m ) is restricted to the gray 
area which is limited from above by the upper bound of 
Eq. 1421 and from below by the lower bound of Eq. H43II . The 
two curves meet for z = k integer. Here m — 2 and n = 0.3. 

The Renyi entropy of Eq. (0J is decreasing function 
of the parameter z. As a matter of fact, it obeys the 
inequality [13J 



I z' 

S z (p) > - 



-S Z '{p) , 



(40) 



which applies for any z ^ z' and for all p. This property 
and the monotonicity of S z (p) respect to the norm of 
Eq. can be used to derive the relation \\p\\ z ^ ||p||z' 
which, when applied to the output state of a channel M., 
implies 



v z {M® m ) < v z ,(M® m ) 



z > z' > 1 



(41) 



In the case of the channel Af n , by choosing z' = k integer 
and using the identity we obtain the upper bound 
for all z > k, i.e. 



1 



(n + l) k — n k 



m/k 



(42) 



In this section we show that the results obtained for the 
channel M n apply also to other Gaussian channel mod- 
els. The first group we analyze is formed by the classical 
channels Q where, as in the case of M n the photonic sig- 
nal from the sender is displaced randomly in phase space 
according to a Gaussian distribution. The second group 
is formed by the channels C where the input signal is lin- 
early coupled to an external environment prepared in a 
Gaussian state. 



A. Classical channels 

Consider CP map Q which transforms the input state 
p into the output state 



Q{p) = / d\ 



„ exp^c-r-C 1 ) 



7r/(2ydrt[T]) 
x exp [(a, a') ■ C^] P exp [— (a, a') ■ C^] 

where £ = (/x, — p*) and 



(44) 



u v 

V u 



u ^ \v\ 



(45) 



is a 2 x 2 positive Hermitian matrix [l^. For T = l/(2n) 
the map Q gives M n of Sec. [I] while for generic Y the 
channel (|44|) is the generalization of Af n to the case of 
non circularly symmetric distribution Q. As shown in 
App. El the map Q{p) can be decomposed according to 
the relation 



G(p) = Et(0 AA n (S(C) p Et(0) E(0 , 



where 



= 1/(2^-1^2) 



-arctanh 



1/2 



+ yju 2 - \v\ 2 



and where 



E(0 ee exp [(r a 2 - i Wf) /2] 



(46) 

(47) 
(48) 

(49) 



is the squeezing operator. In other words, for any input 
state p, the output state Q(p) can be obtained by ap- 
plying the squeezing operator E(£) to p, then sending it 
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through the channel J\f n , and, finally, applying the anti- 
squeezing transformation We can thus consider 

M n as a simplified version of Q where all the squeezing 
operations have been removed. 

An important consequence of Eq. (|46|l is that the z- 
norms of the channels Q and J\f n are identical. In fact, 
using the invariance of the norm under the unitary 
operation £^(£)® m , we can write the z-norm of m uses 
of the channel Q as 

v z {Q® m ) = max ||^ m (£(0® m p£+(£)® m )|U 

p£H® m 

= max \\Ar® m (p)\\ z = p z (Af® m ) (50) 

where, in the second identity, the unitary operator 
£(£)® m has being incorporated in the definition of the 
input state p of the m uses of the channel. In particular, 
for m = 1 and z — k integer, Eqs. (|50() and (|llfl give the 
fc-norm for the single channel use of Q, i.e. 



(2^ 



\2\k 



-i 1/fe 



(1 + 2y/v? - \u\ 2 ) k - 1 



(51) 



According to the decomposition rule of Eq. 146(1 . such 
a maximum is achieved for the anti-squcezed coherent 
states 



(52) 



In fact, feeding the channel Q with this input is equiv- 
alent (apart from an irrelevant unitary transformation) 
to feeding M n with the coherent state |a). Moreover, for 
generic m Eq. 1(50(1 gives 



(1 + 2y/v? - \u\ 2 ) k - 1 



i/k 



(53) 



which proves the Amosov-Holevo- Werner conjecture for 
the channel G, at least for integer z = k. As in the case 
of Eq. 1(52(1 . the input states that achieve the maximum 
((53(1 can be obtained by anti-squeezing the states which 
achieve the maximal output z-norm for the channel A/" n , 
i.e. |«i; ® • • • ® \a m ; -£)m- 



B. Linear-coupling channels 

The linear-coupling channel model represents a com- 
munication line where the input photons (described by 
the annihilation operator a) interact with an external 
environment (with annihilation operator b) through the 
beam splitter unitary operator 



U = exp 



(a) b — ab^) arctan 



1-7? 



V 



which transforms the fields according to 

a — > U^aU 
b — ► U*bU 



y/rja + sj\ - 7] b 



b - \j\-r\ a 



(54) 



(55) 



with 7] e [0, 1] being the beam splitter transmissivity. For 
77 = 1, U is the identity and the input signal is decoupled 
from the environment; for rj — 0, instead, U is a swap 
operator which replaces the input signal with the envi- 
ronment input state. The CP map of the linear-coupling 
channel is obtained by coupling the input state of the sig- 
nal p with the input state of the environment Tb through 
U and then by tracing away the mode b. The resulting 
output state is then 



£{p) = T n [Up® Pb C/t] 



(56) 



If pb is a Gaussian state, the CP map C is Gaussian. 
In what follows we will assume pb to be the squeezed 
thermal state 



p b = Ej(0 r b (n) £ b (£) 



(57) 



where £& and Tb(n) are, respectively, the squeezing oper- 
ator and the thermal state 1)57(1 of the b mode. For the 
channel l(56|l a decomposition rule analogous to Eq. 146(1 
applies, namely (see App. IB 2(1 



C(p) = Et(0 £„(£(£) P S f (0) s (0 . 



(58) 



with £ n (p) the CP map 1(56(1 where the environment is in 
the thermal state pb = t(ti). The connection between C 
and £ n is thus analogous to the connection between Q and 
J\f n . In particular we can derive the following identity 



(59) 



which applies for all m integer and for all z. Proving the 
Amosov-Holevo- Werner conjecture for £ n is equivalent to 
proving it for C: moreover, the input states which achieve 
the maximum output z-norm for C are obtained by anti- 
squeezing the input states which achieve the maximum 
for £ n . 

The channel £ n has been extensively studied in 
Ref. 01 where it was shown that it satisfies the relation 

£n(p) = (-A/(i-77)ti £0) (p) = J^(i-n)n ( £o(p)) , (60) 

with £q being the lossy map, where the input photons in- 
teract with the vacuum state of the environment. Equa- 
tion ((60(1 shows that the output of the channel £ n can be 
obtained first applying the lossy map to the input state 
p and then feeding it into the classical channel Af n ■ This 
composition rule has two important consequences. On 
one hand, it implies 

M^r ) - MMi-r,)* ° torn < *, (^a n)n ) ■ ( 61 ) 

In fact, the maximization implicit in the second term is 
performed on a set of input states which form a proper 
subset of the input states which enter in the maximiza- 
tion of the third term. On the other hand, since the lossy 
channel maps coherent input states into coherent outputs 
according to the transformation [3 



£ (\a)(a\) = \y/rja)(y/rja\ , 



(62) 
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equations l|60|) and (|10f> show that when £ n and Af(i- V )n 
act on coherent inputs they produce the same output 
z-norm. This is sufficient to prove that, at least for z = k 
integer, the inequality in Eq. (|61[) is replaced by an iden- 
tity: we have already established in fact that the maxi- 
mum fc-norm of the channel Af n is achieved for a coherent 
state. Hence we can establish the following identity 



where we have introduced the complex linear vector 



1 



[(1 - rj)n+ l] fe - [(1 - r])n} k 



m/k 



,(63) 



which, analogously to Eq. shows that the Amosov- 
Holevo- Werner conjecture applies for the channel £ n at 
least for all integer k. Moreover we know that, as in 
the case of Af n , tensor products of coherent states are 
sufficient to achieve the maximum of Eq. (|63|) . 



III. CONCLUSION 

In this paper we have studied various model of Gaus- 
sian bosonic channel (i.e. the classical maps Q of Sec. 
Ill Al and the linear coupling maps C of Sec. Ill B[) and we 
have shown that the Amosov-Holevo- Werner conjecture 
@ applies to them at least in the case of z — k inte- 
ger. In particular we have proved that tensor products 
of squeezed coherent states are the inputs that achieve 
the maximum output /c-norm for the m successive uses 
of these channels. In the case of the circularly symmet- 
ric channels Af n and £ n the optimal state are just tensor 
product of coherent states. These properties imply that, 
for all integer order greater than 2, the Renyi entropies 
at the output of the channels Q and C are additive, and 
suggest that the same behavior should apply also to all 
the other orders (see Sec. II C|l . In particular, it seems 
reasonable to believe that these channel posses an addi- 
tive Holevo information [l5ll20| . 



APPENDIX A: PROPERTIES OF THE 
OPERATOR 6. 

In this appendix we prove that the operator of 
Eq. (|22[) has the tensor product structure of Eq. i|23|) 
and we derive the identity 



1. Derivation of Eq. lt23ll. 

By using the property 

Dl{n)D s (v)=D s {v-fl) exp[(^. A /t_ A 2. i 7t)/ 2 ](Al) 

of the multi-mode displacement operator defined in 
Eq. (jT%|) . the expression of 6 yields 
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d g e -/?-C-M t+^.Gt-a t-a-G'V f /^j 



A* = (fill $21 

= (Mil)--- 



, ^2ml Mfcl, 



which has km elements, and 

a = (fill a 2 ; • • • ; a k ) 
= (an, • • • , o-iml 0,21, ' 



, 0,2ml a kl, 



(A3) 



(A4) 

5 &fc m ) 



where a sr is the annihilation operator associated with the 
sth copy of the rth use of the channel. In Eq. (| A2|) . G 
and 



i 

n 



(A5) 



are now mk x mk real matrices (1 is the mk x mk iden- 
tity), which are obtained, respectively, by tensoring to 
the mth power the matrices G and C of Eqs. Ij27(l and 
(|26|l . In particular for k ^ 3, G and A have the block 
form 



G' y 



— j^iim 



" -1 


1 


■• 











-1 


1 •• 














-1 •• 














•• 


• -1 


11 


1 





•• 





-1 


■ 


-1 


•• 


■ 


1 ■ 


1 





-1 •■ 


• 








11 


•• 


• 











•• 


• - 


-1 


. -1 





•• 


• 1 


0. 



(A6) 



(A7) 



where now 1 and are the mxm identity and null matrix 
respectively. On one hand these equations show that the 
Gaussian in the integral (|A2|) couples together all the 
annihilation operators a sr which have the same index r 
(i.e. the operators a\ r , a,2 r , • ■ ■ and a,k, r which enter in 
the definition of the vector a r of Eq. I|24|ll. On the other 
hand, Eqs. (|A6|) and (|A7() show that any two annihilation 
operators a sr and a S ' r i with r ^ r' are not coupled by the 
integral (|A2|) . We can hence write such an expression as 
a (tensor) product of m independent Gaussian integrals 
where only the modes associated with a r enters: by doing 
this and using the cyclic symmetry of the matrices G and 
A wc finally obtain Eq. (J23J. 



2. Derivation of Eq. (l33l). 

The identity has been derived in by showing 
that the operator on the left-hand-side of this expression 
has the same symmetric characteristic function of the op- 
erator Qj r of Eq. (|32|l . Here, instead, we show that these 
operator coincide by calculating their matrix elements 
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in the Fock basis of the annihilation operator bj r . For 
the sake of simplicity in the following we will omit the 
indexes j and r. 

Given \p) and \q) Fock states of b with p ^ q, consider 
the quantity 



where u and v are the elements of T defined in Eq. J7J • 
The matrix B has determinant equal to 1 and inverse 



B~ 



a 



-(3* 



a 



(B3) 



( P \e\ q ) 



d"v- 



irniei 



(p\D(— v)\q) . (A8) which diagonalizes T through the relation 



Following the derivation given in p"?| the matrix element 
in the integral can be expressed in terms of the Laguerre 
polynomials as 



(p\D(-v)\q) = 



(A9) 



Replacing this expression in fKSl and using the identi- 
ties m 



dx e-i x L^(x) = ^ 7 Rc 7 > , (A10) 



and 



(AH) 



where 5 pq is the Kronecker delta, we finally obtain 



(p\0\q) 



2/n 



2d 



2d+ \e\ 2 \2d+ |e| 2 
which proves the thesis. 



(A12) 



APPENDIX B: DECOMPOSITION OF THE 
MAPS Q AND £ 

In this section we derive the decomposition rule of 
Eqs. I|46(l and l|58|) which allows to express the map Q 
of Eq. l|4"4")l in terms of Af n and the map C of Eq. (|5f))) in 
terms of £ n , respectively. 



1. Derivation of Eq. (f4?Tl) 

Consider the Hermitian matrix 

B = 



a (3* 
P a 



(Bl) 



with 



l+y/u 2 - H 2 

2^/u 2 - \v\ 2 ' 



1/2 



y/u' 2 - Rj 

\v\ V 2y/u 2 - \v\ 2 



f> = A ' " 



1/2 



B- 1 T ■ B~ 



2 _ U,|2 





\/ u 2 — |u| 2 



.(B4) 



Moreover, when applied to (a, a') this matrix produces 
the Bogoliubov transformation 



= Et(0(o,ot)E(0, 



(B5) 



where E(£) is the squeezing operator defined in Eq. I|49|l . 
Using these properties we can obtain Eq. I| Kill from 
Eq. (|44|1 by performing a change of integration variables. 
In fact, for £ — > £ • i? we have 



s(/0 = / <fc 



j2 A exp[-C • (i?- 1 • r • b- 1 ) ■ ct] 



V(2 v / det[f]) 
exp^cct)-^] p expf-^ct)-^] (B6) 



= / d 2 a 



exp[— 2yf% 



M 2 H 21 



7r(2V« 2 -H 2 ) 
x Et(OZ>(//)E(OpEt(OZ> t (M)E(0, 

which, according to Eq. (JJJ and (|47|l coincides with the 
left-hand-side of Eq. 



2. Derivation of Eq. (|58|l. 

For the sake of clarity, in what follows the operators 
which acts only on the environment will have the sub- 
script b while the operators which act only on the input 
state will have the subscript a. Using the relations (|55j) it 
is easy to show that the coupling operator U transforms 
a 2 + b 2 into itself, i.e. that it commutes with the op- 
erator E (£)E{,(£) which squeezes both the signal mode 
a and the environment mode b by the same quantity 
£. Inserting the identity decomposition Et(£)E a (£) = 1 
in Eq. (|56[) and using the invariance of the trace under 
cyclic permutation, the above property allows us to write 
Eq. 123) as 



= Et(OTr 6 [C/ (£a(0p£t(0®7i(n)) U ^ 

= Et (£) £„(E a (£)pEt (£)) Eo (0 , (B7) 



(B2) which proves the thesis JSHJ- 



9 



[1] P. W. Shor , eprint |quant -ph/0305035 A. S. Holevo, 

eprint quant-ph/0306196 
[2] M. A. Nielsen and I. L. Chuang, Quantum Computation 

and Quantum Information (Cambridge University Press, 

Cambridge, 2000). 
[3] C. H. Bennett and P. W. Shor, IEEE Trans. Inf. Theory 

44, 2724 (1998). 
[4] A. S. Holevo, Probl. Inf. Trans. 9, 177 (1973). 
[5] A. S. Holevo, IEEE Trans. Inf. Theory 44, 269 (1998); P. 

Hausladen, R. Jozsa, B. Schumacher, M. Westmoreland 

and W. K. Wootters, Phys. Rev. A 54, 1869 (1996); B. 

Schumacher and M. D. Westmoreland, Phys. Rev. A 56, 

131 (1989). 

[6] C. H. Bennett, D. P. Di Vincenzo, and J. A. Smolin, Phys. 

Rev. Lett. 78, 3217 (1997). 
[7] V. Giovannetti, S. Guha, S. Lloyd, L. Maccone, J. H. 

Shapiro, H. P. Yuen, Phys. Rev. Lett. 92, 027902 (2004). 
[8] C. King, IEEE Trans. Info. Theory, 49, 221 229, (2003). 
[9] C. King and M. B. Ruskai, IEEE Trans. Info. Theory 47, 

192 (2001). 

[10] V. Giovannetti, S. Guha, S. Lloyd, L. Maccone, and J. 

H. Shapiro, Phys. Rev. A accepted for publication, eprint 

quant-ph/0404005 
[11] V. Giovannetti, S. Lloyd, L. Maccone, J. H. Shapiro, and 

B. J. Yen, Phys. Rev. A accepted for publication, eprint 
quant-ph/0404037 

[12] G. G. Amosov, A. S. Holevo, and R. F. Werner, Problems 
Inform. Trans. 36, 305 (2000), eprint math-ph/0003002 
[13] K. Zyczkowski, Open Syst. and Inf. Dyn. 10, 297 (2003); 

C. Beck and F. Schlogl, Thermodynamics of chaotic sys- 



tems (Cambridge University Press, Cambridge, 1993). 
[14] G. G. Amosov and A. S. Holevo, Th. Prob. Appl. 47, 123 
(2001). 

[15] A. S. Holevo, M. Sohma, and O. Hirota, Phys. Rev. A 
59, 1820 (1999); M. Sohma and O. Hirota, Recent Res. 
Devel. Optics, 1, 146-159 (2000) edited by Research Sign- 
post, eprint quant-ph/0 105042 A. S. Holevo and R. F. 
Werner, Phys. Rev. A 63, 032312 (2001); M. Sohma and 
O. Hirota, Phys. Rev. A 65, 022319 (2002). 

[16] U. Grenander and G. Szego, Toeplitz forms and their ap- 
plications, (Univ. of Calif. Press, Berkley, 1958). 

[17] C. Caves, "Hidden- variable model for 

continuous- variable teleportation", in 

http://mfo.phys.unm.edu~caves/reports/cvteleportation. 



[18] D. F. Walls and G. J. Milburn, Quantum Optics (Springer 

Verlag, Berlin, 1994). 
[19] In Eq. H46H the differential form d 2 ( stands for 

d 2 fx = dRe(n) dlm(fi). 
[20] V. Giovannetti, S. Lloyd, L. Maccone, and J. H. Shapiro, 

"Classical capacity of bosonic channels" , unpublished. 
[21] I. S. Gradshteyn, I. M. Ryzhik , Table of integrals, series, 

and products (Academic Press, San Diego, 2000), Ch. 7. 

416. 

[22] The operator O has an orthonormal basis of eigenvectors 
(see for instance Eqs. I23|l . lo 1 1 and i|33pL In general, 
however, its eigenvalues are complex numbers: is thus 
not Hermitian as erroneously reported in the published 
version of this paper (Phys. Rev. A 69, 062307 (2004)). 



